The Gaussian free field in interlacing particle 

systems 



We show that if an interlacing particle system in a two-dimensional 
lattice is a determinantal point process, and the correlation kernel can be 
expressed as a double integral with certain technical assumptions, then 
the moments of the fluctuations of the height function converge to that 
of the Gaussian free field. In particular, we examine a specific particle 
system that was previously studied in and compute the corresponding 
Green's function. 

1 Introduction 

The Gaussian free field is widely considered to be a universal object describing 
the flucutations of heights of random surfaces. Previous work has rigorously 
shown this to be the case in specific models {[l],[in])- In this paper, we show 
that if an interlacing particle system in two dimensions can be described as a 
determinantal point process, and mild assumptions are made about the corre- 
lation kernel, then the covariances of the fluctuations of the height function are 
governed by a particular Green's function. A general formula for the Green's 
function is given. 

In particular, we will use this general theorem to show that in an interlacing 
particle system that arises from the representation theory of the orthogonal 
groups, the Green's function is 



Note that G is the Green's function for the Laplace operator with Dirichlct 
boundary conditions on the set H — D := {z S C : 3z > and \z\ > 1}. It 
turns out that there is a map from the surface to H — D. We will show that 
the fluctuations of the height function converge to a Gaussian process whose 
covariance is given by the puUback of Q under D,. 

Particle System. Now let us describe this particle system, which was the 
initial motivation for this paper. Introduce coordinates on the plane as shown 
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in Figure [T] Denote the horizontal coordinates of all particles with vertical 
coordinate m hy yY^ > > ■ ■ ■ > y]^ , where k = [{m + 1) /2J . There is a wall 
on the left side, which forces > for to odd and > 1 for to even. The 
particles must also satisfy the interlacing conditions y^^i < < 2/™^^ for all 
meaningful values of k and to. 




By visually observing Figure [T] one can see that the particle system can be 
interpreted as a stepped surface. We thus define the height function at a point 
to be the number of particles to the right of that point. 

Define a continuous time Markov chain as follows. The initial condition is 
a single particle configuration where all the particles are as much to the left as 
possible, i.e. = m — 2A; + 1 for all k, m. This is illustrated in the left-most 
iamge in Figure |2] Now let us describe the evolution. We say that a particle 

is blocked on the right if + 1 — y^k~i^ ^^'^ blocked on the left if 

— 1 = y^~^ (if the corresponding particle y^Zi or y™~"^ does not exist, then 
y^ is not blocked). 

Each particle has two exponential clocks of rate all clocks are independent. 
One clock is responsible for the right jumps, while the other is responsible for 
the left jumps. When the clock rings, the particle tries to jump by 1 in the 
corresponding direction. If the particle is blocked, then it stays still. If the 
particle is against the wall (i.e. yVL+i, = 0) and the left jump clock rings, the 

V 2 J 

particle is reflected, and it tries to jump to the right instead. 

When y™ tries to jump to the right (and is not blocked on the right), we 
find the largest r e Z>o U {+00} such that y™^' — y™ + i for < i < r, and 
the jump consists of all particles {y™^*}^_Q moving to the right by 1. Similarly, 
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when tries to jump to the left (and is not blocked on the left), we find the 
largest I € Z>o U {+00} such that y™^j = — j for < j < Z, and the jump 

consists of all particles {v^+jY^^Q moving to the left by 1. 

In other words, the particles with smaller upper indices can be thought of 
as heavier than those with larger upper indices, and the heavier particles block 
and push the lighter ones so that the interlacing conditions are preserved. 



Figure 2: First three jumps 




Figure [2] depicts three possible first jumps: Left clock of y\ rings first (it gets 
reflected by the wall), then right clock of rings, and then left clock of y\ 
again. 

In terms of the underlying stepped surface, the evolution can be described 
by saying that we add possible "sticks" with base 1x1 and arbitrary length of 
a fixed orientation with rate 1/2, remove possible "sticks" with base 1x1 and 
a different orientation with rate 1/2, and the rate of removing sticks that touch 
the left border is doubledEI 

A computer simulation of this dynamics can be found at 
|http : //www.math.caltech.edu/papers/Orth_Plaxich.htinl| 

This particle system falls in the universality class of the Anisotropic Kardar- 
Parisi-Zhang (AKPZ) equation with a wall. The KPZ equation was first intro- 
duced in [8] and is of interest to physicists, see [5]. Similar Markov chains have 
been previously studied in without the wall, and in [TS' with a different 
( "symplectic" ) interaction with the wall. 

More general particle systems. By following the proof, the author real- 
ized that a more general statement could be proved. If a point process in a 
two-dimensional lattice is determinantal, and the correlation kernel can be ex- 
pressed as a double integral with certain technical assumptions (see Definition 
|2.1| below), then the moments of the fluctuations of the height function can be 
governed by a Green's function. The exact statement is Theorem |2.3[ We then 

^This phrase is based on the convention that vT is a figure of a 1 X 1 X 1 cube. If one 
uses the dual convention that this is a cube-shaped hole then the orientations of the sticks to 
be added and removed have to be interchanged, and the tiling representations of the sticks 
change as well. 
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can use this theorem to determine the Green's function for the specific point 
process described above. 

Motivations. There are three reasons for proving the results in the more 
general case. The first reason is that the proofs are not much more difficult. 
The second reason is that it is easier to check the conditions for the general case 
than it is to repeat the full calculations. The third reason is that the general 
result tells us that the formula for the Green's function only depends on G^, and 
therefore should only depend on the horizontal movement of the particles. For 
instance, in the model of 2 , the expression for Q was explicitly computed. In 
the model of [6J, the particles had two different jump rates, and the expression 
for Q was the same. Our general results should imply, for example, that Q will 
be the same for any number of different jump rates. In order for a model to have 
a different expression for there should be some change along the horizontal 
direction. For example, in the model considered in this paper, a reflecting wall 
is added to the model of [2 , and we see that Q is changed. 

The correlation kernels in previous problems ([4], [6]) should satisfy the con- 
ditions of the more general case, although the author has not rigorously checked 
the details. 

Conjectures. There are several statements which should be true, but are 
not pursued in this paper. One of them concerns the one-point fluctuations. 
Namely, after a scaling of \/\og N, the fluctuations of the height function at 
a single point should converge in moments to a standard Gaussian. The proof 
would consist mainly in proving that the variance is proportional to log N, which 
by [12j immediately implies the convergence to a Gaussian. 

The appropriate scaling was first demonstrated in |14j using a renormaliza- 
tion group argument. Using this convergence, one can then show that the con- 
vergence to the Gaussian free field occurs in distribution, not just in moments. 
Analogous (rigorous) statements are found as Theorem 1.2 and Proposition 5.6 

It should also be possible to prove similar results when ti ^ T2, but this was 
also not pursued. 

Organization of Paper. In section [2?T] the theorem for more general particle 



systems (Theorem 2.3 1 is stated. In sections 2.2 and 2.3 the algebraic and 
analytic steps of the proof are given. In section |3j it is shown that the particle 



system described above satisfies the conditions of Theorem 2.3 Section 4.2 deals 
with the necessary asymptotic analysis. 

Acknowledgments. The author would like to thank Alexei Borodin and Mau- 
rice Duits for insightful comments. The author was supported by the National 
Science Foundation's Graduate Research Fellowship Program. 
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2 General Results 



2.1 Statement of the Main Theorem 

Suppose we have a family of point processes on X ZxZ>i which runs over time 
t G [0,00). (Note that these are different co-ordinates from the introduction). 
In other words, at any time t, the system selects a random subset X C X. If 
{x,n) e X, then we say that there is a particle at {x,n). For any k > 1 and 
t > 0, let pI : X'' -> [0, 1] be defined by 

pl{xi,ni, . . .,Xk,nk) 

= P(There is a particle at (xj, rij) at time t for each j = 1, . . . , fc). 

Assume that there is a map if on X x X x [0, 00) such that 

pl{xi,ni, . . .,Xk,nk) dct[K{xi,ni,Xj,nj,t)]i<i^.j<k- (1) 

The maps p^ and K are called the kth correlation function and the correlation 
kernel, respectively. 

A function c on X x X is called a conjugating factor if there exists another 
function C on X such that 

/ r /\ C{x,n) 

c[x,n,x ,n = -— — -. 

C(x, n') 

Note that if c is a conjugating factor, then 

det[K{xi,ni,Xj,nj,t)]i<ij<k = det[c{xi,ni,Xj,nj)K{xi,ni,Xj,nj,t)]i<ij<k- 

(2) 

Two kernels K and K are called conjugate ii K = cK for some conjugating 
factor c. 

If a correlation kernel exists, the point process is called determinantal. On a 
discrete space, a point process is determined uniquely by its correlation functions 
(see e.g. Therefore, if we are given two determinantal point process on a 

discrete space with conjugate kernels, they must have the same law. 

The set Z x {n] is called the nth level. Given a subset X C X, let m„ be 
the cardinality of the set X n (Z x {n}). Assume that the numbers m„ take 
constant finite values which are independent of the time parameter t. In words, 
this means that the number of particles on the nth level is always m„. Further 
assume that to„ < m„_|_i < to„ + 1 for all n. Let xj""* > x'^^^ > . . . > Xml 
denote the elements of X D {Z x {n}) . A subset X is called interlacing if 

^k+i^^ < ^fe"'' - ^k'^^^^ when run+i = m„, 

(n+l) ^ (n) ^ (n+1) , , , 

Xk+i <Xk <xl ' , when to„+i = m„ + 1. 

Assume that at any time i, the system almost surely selects an interlacing subset. 
Let 5n equal run+i — rrin- 
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Define the random height function by 

h:Xx R>o Z>o, 

h{x,n,t) = \{{s,n) G X : s > x}\. 

In words, h counts the number of particles to the right of {x,n) at time t. 

We wish to study the large-time asymptotics of this particle system. Let x = 
[iVi/], n = [Nrjjjt = Nt, where iV is a large parameter. Define V cRx M+ x M+ 
to be 

V := {{i/,r],T) : lim p{{x,n) > 0}. 

Let Hn be defined by 

i/jv : K X M+ X M+ ^ M, 

_ffjv(i^, 77, r) :— h{x, n, t) — E/i(x, n, t). 

In words, i?jv is the fluctuation of the height function around its expectation. 

Before stating the theorem, we need to state some more assumptions on the 
kernel. 

Suppose the kernel K is conjugate to a kernel K such that K satisfies the 
following property: There is a number L such that whenever x, x' > L, 

K{x, n, x' , n , t) + K{x, n, x' — 1 + n' + l,t) + K{x, n, x + (5„, n' + 1, i) 

(x,n) = (x\v!) 
lo, otherwise, 

K{x,n, x',n',t) + K{x + 1 — ^„,n — 1, x' , n' , t) + K{x — 6n,n— 1, x', n',t) 

j 1, {x,n) = {x',n') 



0, otherwise. 



(4) 



Further suppose that for x',x" > L, 

K{x, n, x', n', t)K{x", n" ,x — 1 + 5n,n + l,t) ^ as x ^ 00 (5) 

K{x,n,x' ,n' ,t)K{x" ,n" ,x + 5n,n+ — )• as a; ^ 00 (6) 

K{x, n, a; — 1 + n + 1, ^ as a; — )• 00 (7) 

K{x, n,x + 5n,n + l,t) ^ 1 as X ^ 00. (8) 

Suppose G{v,'q,T,z) is a complex-valued function on R x R_|_ x M+ x C. 
To save space, we will sometimes write G{z). Expressions such as G',G,^,G'^ 
will be shorthand for dG/dz, dG/dv and d^G/dzdv, respectively. Assume 
= G(2;).Also suppose there exists a differentiable map f2 from T) to the 
upper half-plane H = {2 e C : '^{z) > 0} such that f2 is a critical point of G. 
In other words, 

G'{u, ri, T, n{u, 7], r)) = for all {u, rj, r) e V. (9) 
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Note that 17 need not be onto. For any (77, r), if the set {i^ G M : {v, rj, t) G D} 
is nonempty, let 92 ('y,''') denote its supremum. 

Definition 2.1. With the notation above, a determinant al point process on 
Z X Z>i is normal if all of the following hold: 

• For all yy, r > 0, the limit fl{q2{ri, t) — 0, 77, t) exists and is a positive real 
number. 

• For all 77. r > 0, as approaches g2('7, t) from the left, G"{i', rj, r, il{L', ry, r)) - 
0(^(77, r)-;/)i/2). 

• K is conjugate to some K such that ([3|-([8| hold for some integer L. 

• Set t = Nt, Xj = [Nvj] and nj — [Nrjj] for j ~ 1, 2, where {i^j, rij,T) £ T). 
Let U,j denote Q{i>j,rij,T) and let Gj{z) denote G{i'j,r]j,T, z). For finite 
integers ki and fc2, as — ^ c», 



l\ f f exp(iVG(i/i,77i,r,z)) / , , „ , , , ^o/„JV« 



i4r(xi, 7ii + ki,X2,n2 + k2,t) 

Stti/ Jra exp(iVG(7^2,??2,-r,u))) 

where Fi and r2 are steepest descent paths, k < ^{Gi{fli) — 02(^2)), and 
fkmkn are complex- valued meromorphic functions satisfying the identity 

fkik2{zi,Z2)fk2k3{z2,Z3) . ■ . fk,^_,kA^r~l, Zr)fkr.kdZr, Zl) 

= f{zi,Z2)f{z2, Z3)... f(Zr-l,Zr)f{Zr,Zi). 

Here, we have written / for /qq. 

• For any Z > 3, the following indefinite integral satisfies 

^j^ /(M;),M;|i)) ,,^,0, (11) 

where the sum is taken over all Z-cycles in Si and the indices are taken 
cyclically. 

• For any finite interval [a,b], G E C'^[a,b] and the Lesbesgue measure of 
the set {x G [a, b] : I'{x) G 27rZ + [—5, 6]] is ©((J") for some positive a. 

The following remarks will help explain the definition. 

Remark 2.2. (1) The assumption that Q,{q2{ri, r) — 0, 77, r) > occurs naturally. 
One often finds that for fci 7^ fc2 G Z, 

K{[Nv]+ki,[NT^],[Nv]+k2,[Nri\,NT)- ' - ^ 



N^oo J '^'^ '^'^ ' 2TriJa zfei-'^s+i 7r(fc2-fci) 
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where the contour crosses the positive real hne. By setting = e'f , we see that 
the right hand side reduces to the ubiquitous sine kernel. When ki ~ k2 = 0, 
we see that 

lim prUNi^i [Nv]) = Uiogn- logn) = "'"g"^"'^'"^ 

N^oa Zn TT 

Since the left hand side equals zero, we expect arg fl(q2{r], r) — 0, ij, r) = 0. 

(2) Since G{z) = G{z), this means that G has a critical point at both and 
fj. As r2 approaches the real line, the two critical points coalesce into a triple 
zero, so G"{t,ri,T) converges to as t approaches q2{r],T). We need a control 
for how quickly this convergence to occurs, in order to order to control the 
behavior near the boundary of V. More specifically, it controls the bound in 
Proposition |4.f 1| 

There is a heuristic understanding for why (2) should hold. The function G 
has two critical points which coalesce into a triple zero. The simplest example 
of such a function is G{t,x) = — tx as t approaches 0. In this case, the 
solution to G'{t,x) = is VL{t,x) = t^/"^. Then G"{t,n{t,x)) = 2t^/^. 

(3) Assumptions ([3])-([8]) will be elucidated when we interpret the particles 



as lozenges. In particular, see remark 2.6 



(4) It is common for the kernel to be expressed in this form; previous ex- 
amples are |1] and [5]. If the kernel has a different expression with the same 
asymptotics as in Propositions |4.5| and |4. 11[ the results still hold. 



(5) In particular, (111 holds if there always exist M-substitutions and an 
expression Y such that 

where = zi and Wi+i = ui. This is because of Lemma 7.3 in [7], which 
refers back to [T] , which says that 

' 1 



(6) This is a technical lemma which allows Lemma 4.1 to be applied. 
We can now state the main theorem. 

Theorem 2.3. Suppose we are given a normal determinantal point process. For 
1 ^ J ^ fc; let = (i^ji %j t) be distinct points in V , and let flj = Qii'j, rjj, r). 
Define the function Q on the upper half-plane to he 

n \ f r r f{zi,z2)f{z2,z,) ^ ^ 



The 



lim EiH^i^i) ■■■Hm{xu))-- 

N—t-OO I ^^-^k J — ^ 



k/2 



IS even 



0, k is odd, 
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Figure 3: In this example, the integers m„ equal 1, 1. 1, 2, 3, 3, 4 . . .. The black 
line on the left represents the points where x ~ 0. Examples of x^,"-* are x^f^ = 




x=0 



where 'is the set of all involutions in Sk without fixed points. 

Remark 2.4. We note that these are the moments of a linear family of Gaussian 
random variables: see Appendix A. Using the results of [12], it should be possible 
to show that iJjv (>«■) / \J^8j:Hn [x) converges to a Gaussian, but this was not 
pursued. 



2.2 Algebraic steps in proof of Theorem 2.3 



The most natural way to view X is as a square lattice. However, it turns out 
that a hexagonal lattice is more useful. To obtain the hexagonal lattice, take 
the nth level and shift it to the right by (n + l)/2 — m„. See Figure [sj 

Figure [3] also shows that the particle system can be interpreted as lozenges. 
Each lozenge is a pair of adjacent equilateral triangles. See Figure |4] 

By setting the location of each triangle to be the midpoint of its horizontal 
side, each lozenge can be viewed as a pair (a;, n, a;', n'), where the black triangle 
is located at {x,n) and the white triangle is located at {x' ,n'). For example, in 
Figure|3]there are lozenges (1, 3, 1, 3), (2, 3, 2, 4), and (0, 3, 1, 4). The three types 
of lozenges can be described as follows. For lozenges of type I, (x', n') = {x, n). 
For lozenges of type II, (a;', n') = (x — 1 + (5„, n + 1). For lozenges of type III, 
(x', n') = (a:; + n + 1). Note that a lozenge of type I is just a particle. 

We say that (x, n, x' , n') e X x X is viable if (a;', n') = (a;, n), (x — 1 + 5„, n + 
1), or (x + (5„,n + l). A sequence (xi, ni, x'j^, n'j^), . . . , (xfe, n/j, x'^, n'^) of viable 
elements is non- overlapping if (xi,ni), . . . , (xfe,rtfe) are all distinct from each 
other and (x'^, n'^ ),..., (xj,, n'^.) are also all distinct from each other. We do, 
however, allow the possibility of (xi^rii) — (x^-,n^). 
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Figure 4: Lozenges of types I, II, and III, respectively. Note that lozenges of 
type I occur exactly at the same places as particles. 



^ n 



The statement and proof of the next proposition are similar to Theorem 5.1 
of kn:BF. 

Proposition 2.5. Suppose the kernel K is conjugate to some K such that 
([8| hold for some L. If t > 0,xi,x[, . . . ,Xk,x'i^ > L, and {xi,ni, x'i,n[), . . . , 
(xk,nk,x'i^,n'f,) is a sequence of non- overlapping viable elements of X x X, then 

V{There is a lozenge {xj ) at time t for each j — 1, . . . ,k) 

= det[K{x^,ni,x'j,n'j,t)]i<i^j<k. (12) 

Remark 2.6. The equations ([3|~([8| can now be intuitively understood. Equation 
([3| says that each black triangle is located in exactly one of the three lozenges 
around it, and equation Q makes an identical statement for white triangles. 
Equations ^ and ([7| say that lozenges of type II almost surely do not occur 
far to the right of the particles, with ^ controlling the off-diagonal entries in 
the determinant and ([7| controlUing the diagonal entries. Similarly, equations 
([6| and ([s]) says that lozenges of type III almost surely do occur far to the right 
of the particles. This intuition will be exploited in the proof of Thereom|2.5| 



Proof. We proceed by induction on the number of lozenges that are not of type 
I. When this number is zero, the statement reduces to ([T]) and ([2|. 

For any set S = {{xi,ni,x'i,n'i), . . . ,{xk,nk,x'j^,n'f.)} of non-overlapping. 



viable elements, let P{S) and D{S) denote the left and right hand sides of (12), 
respectively. First, as a preliminary statement, we prove that if {xk+i,nk+i) ^ 
{Xr,nr) for 1 < r < k, then 

D{SU{{xk+i,nk+i,Xk+i,nk+i)})+D{SU{{xk+i,nk+i,Xk+i-l+Sn, ?^fe+l+l)}) 
+ D{S U {{xk+i,nk+uxk+i + 6n,nk+i + 1)}) = D{S). (13) 



Note that if D is replaced by P in (13), the statement is immediate, since the 



black triangle at {xk+i,nk+i) must be contained in exactly one lozenge. 

Expand the determinants into a sum over permutations a € Sk+i- Now split 
the sum into two parts, one part where the a fixes fc -I- 1 and one part where a 
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does not fix fc + 1 . The first part equals 



a-{k+l)=k+l r=l 

X ^K{xk+i,nk+i,Xk+i,nk+i,t) + K{xk+i,nk+i, Xk+i - 1 + + l,t) 

+ K{xk+i,nk+i,Xk+i + Sn^^^,nk+i + l,t) , 
wliich, by equals D{S). The second part equals 

k 

^ sgn(cr) Jl 

cr(k+l)^k+l r=l 

+ K{x„(k+i),'n^{k+i),Xk+i + Sn^_^i,nk+i + l,t) . 

Since (a;o-(fe+i), "■^(fc+i)) 7^ (a^fc+i, "fc+i) by assumption, ^ implies that this 
expression equals 0. Therefore (13) holds. 

In a similar manner, if {x'^^i, 7^ {x'^, nj.) for 1 < r < k, then Q implies 

that 

L'(S'U{(xfe+i,nfc+i,2:fc+i,nfc+i)})+D(S'U{(xfc+i+l-5„,nfe+i-l,a;fc+i,nfc+i)}) 
+ i?(5U{(xfe+i - (5„,nfe+i - l,a:fe+i,nfe+i)}) = 0(3). (14) 

Again, the statement holds if D is replaced by P. 

Now let us prove that D and P still agree if we add a lozenge of type II to S. 
Suppose that (x, n, x— 1+(5„, rt + 1) is viable and that SU{{x, n, x— 1+(5„, n + 1)} 
is non-overlapping. Then equation (|13|) is equivalent to 



U {(x, n, X - 1 + (5„, n + 1)}) 
= D{S) - D{SU{{x,n,x,n)}) - D(S' U {(x, n, x + (5„, n + 1)}), (15) 

and the same holds for P instead of D. By the induction hypothesis, 

D{S) = P{S), 
D{S[J{{x,n,x,n)}) = ^(5" U {(x, n, x, n)}) 
D{S U {(x + Sn, n + 1, X + 5n,n + 1)}) = P{S (J {{x + 6„,n + l,x + S„,n + 1)}). 



Thus, (15) implies 

D{S U {(x, n, X - 1 + (5„, n + 1)}) - P{S U {(x, n, x - 1 + (5„, n + 1)})) 

= -D{S U {(x, n, X + Sn, n + 1)}) + P{S U {(x, n, x + (5„, n + 1)}). 
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Assume for now that (xj,, n'^) ^ (.x + (5„, n + 1) for 1 < r <k. Then we cam 
apply equation (14 1, which implies that 

D{S U {{x, n, X - 1 + 5„, n + 1)}) ^ D{S) 

- D{S U {{x + Sn, n + 1, x + (5„, n + 1)}) - D{S U {(x + 1, n, x + n + 1)}), 

(16) 

and the same statement holds for P. Thus, 

- D{S U {(x, n,x + 6n,n + 1)}) + P{S U {(x, n, x + (5„, n + 1)}) 
= D{S U {(x + 1, n, X + (5„, n + 1)}) - P{S U {(x + 1, n, x + (5„, n + 1)}). 



If S* U {(x + l,n,x + (5„,n + 1)} is non-overlappinng, then (15 1 is again 
applicable. We repeatedly apply (jlSl) and ( 16 1 as often as possible. First, 



suppose that this can be done indefinitely. Then 



\D{S\J {{x,n,x - 1 + 5„,n + 1)}) - P(S'U {(x,n,x - 1 + (5„,n+ 1)})| 
= lim \D{Syj{{x+M, n, x-l+S^+M, n+l)})-P(5'U{(x+Af, n, x-l+S^+M, n+l)})|. 

A/— voo 

Since lozenges of type II almost surely do not appear when we look far to the 
right of the particles, 

lim F(5U{(x + M,n,x- l + 5„ + Af,n + l)}) =0. 

M— >-oo 

By expanding the determinant into a sum over S'fc+i, ([5| and ([7| imply that 
lim D{S\J{{x + M,n,t,x- 1 + Sn + M,n+ 1)}) = 0. 

M — 7-co 



Now suppose that ( |15[ ) and ( 16 1 can only be applied finitely many times. This 
means that D{S U {(x, n, x - 1 + (5„, n + 1)}) - P{S U {(x, n, x - 1 + (5„, n + 1)}) 
equals either 

£>(S'U {(x + M, n,x + M + 6n, n + 1)}) - P{SU{{x + M, n, x + M + (5„, ?i + 1)}) 



or 



D{SU{{x+M+l,n,x+M+Sn,n+l)})-P{SU{{x+M+l,n,x+M+Sn,n+l)}) 

In the first case, S'U{(x + Af + l, n, x + M + 5n, n+1)} is non non-overlapping. 
This implies D{SU {{x + M + l,n,x + M + Sn,n + 1)}) = (because two of the 
rows are idential) and P{S U {(x + M + l,n,x + M + 6n,n + 1)}) = (because 
a triangle cannot be in two different lozenges at the same time). Thus, D and 
P agree. A similar argument holds in the second case. Thus, D and P agree 
whenever a lozenge of type II is added to S. 

An identical holds for type III lozenges, except that we use ^ and (H 
instead of ^ and ([?]). □ 
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We have been describing a lozenge as a pair {x,n,x' ,n'). It can also be 
described as {x',n',X), where {x',n') is the location of the white triangle and 
A G {I, II, III} is the type of the loznege. Thus the proposition can be restated 
as the following statement. 

Corollary 2.7. For any non-overlappign {x'l, n[, Ai), . . . , (xj,, nj,, Afc), 

V{There is a lozenge [x'j, n'j, Xj) at time t for each j = I, . . . ,k) 

= det[Kxixi,ni,Xi,x'j,n'j,t)]i<^j<k, 

where 

{K{x,n, x' ,n' ,t), when X ~ I 
K{x — 5n-i,n — l,x' ,n' ,t), when X — II 
K{x — 5n-i — 1, n — l,x\n\t), when A = /// 

Proof. This is a result of the correspondences 

{x',n\I) iff {x',n\x\n'), 
(x', n' , II) iff [x' — i5ri'-i, n' — 1, x' , n'), 
(x', n', ///) iff {x' - 5n'-i - 1, n' - 1, x' , n'). 

□ 

There are two different formulas for the height function. One formula is 

h{x, n) — ^ l(lozenge of type I at (s, n)). (17) 



It is possible to only use (17) to complete the proof. However, when there are 
multiple points on one level, i.e. not all rji, ... ,rik are distinct, the computation 
becomes much more complicated. This is because lozenges of type I will appear 
in multiple sums of the form (17). We can avoid this difficulty by introducing 
another formula for the height function: 

h{x,n) = h{x + 5n + (5„+i + . . . + Sn'-i,n') + i?„,„/(a;), (18) 

where, for n < n' , 

n 

IIn,n'{x) = - ^ 1 (lozenge of type II at (x + (5„ + (5„+i ...+ (5p_i,p)). (19) 

p— n+1 

Therefore, the expression 

E n[[hix,,n,)~E{h{x,,n,))]] (20) 
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can be expressed as a sum of terms of the form 

(k' k 
j=i i=k'+i 
Lemma 2.8. Assume that the following sets are disjoint: 

{(s, rij) : s > Xj}, l< i <k' 
{{xi + (5„, + (5„,+i . . . + 5p-i,p) --ni + l <p < n'l}, k' + 1 < I < k. 

Then 



(21) 



(21) 



E ••• E E ••• E det 



Si >xi 

where the matrix blocks are 



All Ai2 
A21 A22 



All = [(1 - 5ij)K{si,n,,Sj,nj,t)]i<ij<k' 

M2 = {K{si,ni,Xri,-pj,t)\i<i<f , k'+i<j<k 

A21 = [-K{xi - l,Sj,nj,t)]k'+i<i<k, i<j<k' 

A22 = [-(1 - Sij)K{xi - Sp^^i,pi - l,Xj,Pj,t)]k'+i<i^j<k 



Proof. By applying Corollary 2.7 to (17) and (19 1, we see that 

(k' k 
'[[h(xj,nj) Yl Hn,^n',{xi) 
i=l l=k'+l 



(22) 



equals the right hand side of (22) with the (1 — 6ij) terms removed. It is well- 
known that subtracting the expectation corresponds to putting zeroes on the 
diagonal. For example, this is noticed in the proof of Theorem 7.2 of [7]. □ 



Write the determinant in (22 1 as a sum over permutations a in S^- If the 



cycle decomposition of a contains the cycle (ci C2 ... c^) of length r and M 
denotes the matrix in the right hand side of (22), then the contribution from a 
is 

E • • • E E • • • E sgn(a)Af,,,,Me,e3 • • • M,^,, {■■■){■■■), 

where (•••)(■■■) correspond to other cycles of cr. Let ipc, denote Sc, if 1 < c,, < 
/c', and if k' < < k. Since the sum over tpc^ only affects the matrix terms 
M=^_ic^ and Mc^c+n the contribution from u is 



Me,,, (...), 



(23) 
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where (. . .) denote other cycles. In other words, the contribution from a can be 
expressed as a product over the cycles in the cycle decomposition of a. 

Note that if a fixes any points, then the correponding contribution is zero 
because all the diagonal entries are zero. 



2.3 Analysis steps in proof of Theorem 2.3 



In (22]), set Xj = = [Nrji], and t — Nr. Our goal is to find the limit of 

(22) as — oo. Expanding the determinant into a sum over a G 5*^, we just 
saw that the contribution from a fixed a is of the form (23). First note that if 
any of the denotes pc- , then 



.Mr 



0. 



This is because each Mc-cj^i is proportional to 1/iV (by Proposition 



4.5 



so 



. Mc,,ci is proportional to N^'"" , but the sum is only taken over 
0{N^~^) terms. Therefore, (20) can be expressed as a single term of the form 
in (21 1, and in this term k' = k. 



ClC2^"C2C3 

rr—1 



Now we will prove (stated as Theorem 2.10 below) that 



■ Mr 



dzi 



dz, 



GUz- 



Once this is proven, ( 11 ) implies that the total contribution from Sk —J^k equals 
zero. When / = 2, then the right hand side is just 0{ili,il2), completing the 
proof of Theorem 2.3 

Recall the definitions of G and fl from section 2.1 Set 6* : 2? — > [0, vr) to be 

^ii^, V,t)^1: arg G"(i/, t], t, r], r)). 



Proposition 2.9. Fori — 1,2,3, let {vi^rji^r) GV, Xi — [Ni^i], rii — [Nrji] and 
t = Nt. For i = 1,3, let Gi{z) denote G{h'i,r]i,T, z), let 9i denote 9{vi,rji,T) 
and let denote n{iyi,r]i,T). Let r+ :— {r2(i/, 772, r) : 1^2 l£ < 12(112, t)} and 
r_=f+. Let Gliz) ^ {d^/dzdv)G{iy2,T2,T,z). Then 



^ K{xi,ni,y,n2,t)K{y,n2,X3,n3,t) 
y>[Ni^2] 

■ Y \ gAf!R((Gi(Oi)-G3(03))) 
O I — 



^-kN ^\G'[{^i)\^\G'i{'^^)\ 7r+ur_ 2^G',^(z) 



dz 



X /(r!i,z)/(z,r!3) 
+ /(r!i,z)/(z,n3) 



giJVG(Gi(Sli))-iei 

gjArs(G3(n3))+ie3 
giJVG(Gi(ni))-iei 



+ /(f]i,z)/(z,f]3) 



-iAr9'(Gi(Oi))+iei 



g-iAfa(G3(n3))-i6(3 



+ /(f^l,z)/(z,n3) 



glAr9'(G3(03)) + ie3 

g-jAr3(Gi(Oi))+iei 
„-iAra(G3(03))-ie3 



(24) 
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Proof. Let G2{z) denote G{[y/N],ri2,T, z), let 02 denote 9{[y/N],r]2,T) and ^2 
denote il{y/N,r]2,T). Fix some /3 € (—1/2,0) and split up the sum into two 
parts: the first part is from [-/Vi^2j to lN{q2 — while the second sum is 

from [N{q2 - N^)\ to [iVqaJ ■ Since there are no particles to the right of Nq2 
in the limit N — oo, the sum from Nq2 to oo can be ignored. It is common to 
refer to the first sum as the bulk and the second sum as the edge. First examine 



the bulk. By Proposition 4.5 



K{xi,ni,y, 712, t)K{y, ^2, X3, 71,3, t) 

„Ar5f((Gi(Oi)-G2(02))) 



,AfJf((G2(n2)-G3(n3))) 



27riVV|G'/(r!i)|V|G'2'(fl2)| 2nN y\a[{n2y\^/\^[{M 



f{ni,n2)fin2,^3) 



(3 



giiV9(Gi(ai))-iei giAr9(G2(n2))-«e2 
giArG(G2(n2))+i«2 e*^9'(G3(n3))+ie3 

h o{G'^{n2y^N-^) + o{G'^{n2)-^ N-^ 



(25) 



where O denotes the other fifteen terms that occur in the sum. First let us 
examine the error term in the bulk. 

each term in the error is bounded by {N^/'^y^N-^ 



2.1 



By (2) of Definition 

and {N^^^)^'^ N^^ , respectively. There are ~ N terms, and since /? > —1/2, we 
must have -2/3 - 3 + 1 < -1 and -7/3/2 - 4 + 1 < -1. Therefore the sum is 
o{l/N). 

Now let us return to the main term in the bulk. For eight of the sixteen terms 
in O, the expression ^^^^(02(^^2)) cancels in the numerator and the denominator. 



By Propositions 4.2 these eight terms are o(l/iV). By 4.4 the other eight terms 
equal 



,AfK((Gi(f2i)-G3(03))) 



-2ie 



2nNy^\GUMV\G[[iM 

X f{ni,n2)f{n2,n3) 



2n\G'^{n2)\ 

giAfS(Gi(ai))-iei 

giJVS(G3(Jl3))+ie3 



dv 



where . . . represent the other seven terms. Of the eight total terms, four have 
/(•, r22)/(^^2, •) E^nd four have /(•, r22)/(^^2, •)■ For the four terms with the 
expression ^2, make the substitution z = ^{v, 772, ■'')• The new integration path 
is T^. By taking the partial of ^ with respect to v and using the chain rule. 



which implies 



-2l02 



an 

dv 



rdv = 



Gim 

G"(17)' 



dv 



dz 



2^|G^'(r!2)| 2^G^'(r!2) 2^GUz)' 
For the four terms with fl2, make the substitution z — Q{v,ri2,T). The path of 
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integration is r_. Finally, the integral becomes 

I \ gAfSR((Gi(Oi)-G3(03))) r 



Z)/(Z, ^3)^,^^^G;^^^^^y^:if^ + f{Cli,z)f{z, »3) g,jV9(G3(03))+.g3 

giJVG(Gi(ai))-i6ii _ _ g-iArQ(Gi(ni))+ie- 



g-iJVQ(G3(n3))-ie3 ' ^ ^ ' ■^^g-j7VS(G3(03))-ie3 



Now we sum over the edge. By Proposition 4.11 and (2) of Definition 2.1 
the sum is bounded above by 

q2N + 1 

y={q2-NP)N y=0 

As long as /3 < 0, the sum over the edge is also o{l/N). □ 

Theorem 2.10. For i — let (vi^rji^T) G T) and set Xi = [Nvi],ni = 

Nrji. Fori = let Gi{z) denote G{i'i,rii, z), let 9i denote 9{vi,rii,T) 

and let $7,; denote ri[ui,rji^T). Let :— {^l{v^rii,T) '■ vi < v < '?2('7i, ''')} (ind 
T- = f+. Then 

I 

Vi>[Nvx\ yi>[Nvi\i=l 



2tx ) Jr+ur- -'r+ur- G'^(zi) Gl{zi) 

The indices are taken cyclically. 



Proof. By Proposition 4.5 the product has 4' terms. Each application of Propo- 
sition [279| decreases the number of terms by a factor of 4, so repeated applications 
of Proposition |2.9| yields the result. □ 



3 Specific Results 

3.1 Particle system with a wall 

We now return to the particle system with a reflecting wall described in the 
Introduction. For notational reasons, it is more convenient to use different co- 
ordinates. Instead of labeling the levels as 1, 2, 3, . . ., it is more convenient to 
label them as (1, -1/2), (1, 1/2), (2, -1/2), (2, 1/2), . . .. If the (ni,ai) is at least 
as high as the (^2,02) level, then this will be denoted as (ni,ai) > (712,02). 
This happens if and only if 2ni + ai > 2n2 + 02. Using the notation of Section 
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2.2 m(^n,a) — " E^nd ^(,i,a) =0 + 1/2. Along the horizontal direction, we will use 
a square lattice, so that the particles live on N instead of 2N or 2N + 1. 
Let niai (dz) be defined by 



dz 
2iz ' 

-{z^/^ ~ z-^/^)^dz 
Mz 



ai=-l/2, 
, ai = l/2. 



Let ji^^^^' denote the (normalized) Jacobi polynomial with parameters 
(±1/2, —1/2). The normalization is set so that for any nonzero complex number 

z, jI*"'"'^^' '^^'^^ satisfies 



,(-1/2,-1/2) 
-'s 

j(l/2,-l/2) 



^5+1/2 _ ^-s-1/2 



,1/2 



^1/2 



(26) 
(27) 



Let W'^°'' ^Z^-* (s) be defined for nonnegative integers s by 

'2, if s > 0, a = - 5, 
^(a,-l/2)(^) ^ ^ 1, if s = 0,a = -i, 

1, if s > 0,a = i. 



Note that for a = ±1/2, 



j(a -1/2) 



ma{dz) = Ss 



(28) 

By Theorem 4.1 of 3], the correlation functions are determinantal with ker- 



nel 



K{ni, ai, Si,n2,a2, S2,t) 



z-\-z 



ly 



(i^+|— !- - 1)"2 z + z-^ - V - V- 



—■ma^{dz)dv 



(29) 



- 1 



ma-,{dz) , (30) 
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where the z-contour is the unit circle and the w-contour is a circle centered at 
the origin with radius bigger than 1. 



Theorem 3.1. The determinantal point process is normal. The Green's func- 
tion is given by 

Giz, w) = ^ log ( ^^^_^ 

Once we prove the point process is normal, the expression for the Green's 
function follows from Theorem 12.31 with 

G[v,r],T;u) = t h^ylog I 1 1 - i/logu, 

f[u,v) = — — 33-. 

In section [3^ we show that the third condition in Definition |2.1| is satisfied. In 
section [331 show that the fourth and second conditions are satisfied. Since 
these are conditions are the hardest to prove, we will focus mainly on their 
proofs. The fifth conditions follows from the substitution Uj = Zj + zj^ and (5) 
of Remark [221 



3.2 Algebraic steps in proof of theorem 3.1 

Proposition 3.2. LetCo{n,a,s) equal 



Co(n, a, s) 



n-l 



(-l)^(-2) 
(-l)^(-2)", a -1/2 



-1/2 



and CQ{ni,ai,si,n2,a2,S2) = Co(ni, ai, si)/Co(n2, 02, S2). Then K = cqK sat- 
isfies for L^l. 



Proof. Using (26)-(27) and the orthogonality relation (28 1, it is straightforward 
to check that ^3^ and ^ hold. What happens is that in the left hand side of 
([3| or Q, one obtains six terms, three of which come from (291 and three of 
which coe from (30). The three terms from (29) always sum to 0, while the 
three terms from (30) sum to or 1. 

Now we will prove ( [7|-([8l ) when ai = —1/2. The term (30) equals zero, so 
we only need to look at ( 29 ) . Explicitly, the expression is 



K{n, -1/2, s,n, 1/2, s',t) 




2iz ' 
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and we want the asymptotic result when s, s' — > oo in such a way that s — s' is 
or 1. Expand the paranthetical expression v'^ — to get two terms, 

each of which is a double integral. Since 1 = |z| < the term with v^^ ^^1"^ 
goes to zero. For the remaining term, expand + z^^ to get two terms. For 
the term with z'*, make the substitution z^^ . What remains is 



e 



27r2j / 7|^|=i gt(^^±^) z^ V - \ ( "+^ ' - 1)" z + z-i - w - v-i 2zz 

Now deform the z-contour to the circle |z| = 1 + 2e and the w-contour to the 
circle \v\ = 1 + e, where e > 0. With these deformations, |?;| < |z|, so the double 
integral goes to zero. However, residues are picked up when the contours pass 
through each other. These residues equal 

„'_„ z dz 
z 

|z| = l+2e Z — 1 2iz 

There is a residue ai z — 1 which equals —2, and a residue at z = which equals 
for s' > s and 2 for s > s' . Since CQ{n, —1/2, s,n, 1/2, s) = —1/2, this proves 
([7| and ([s]) when ai = —1/2. The case when ai — 1/2 is similar. 

It remains to show ^ and Q . When considering the product of two kernels, 
we obtain a quadruple integral. After the substitutions zi i— z^^ and V2 i— > 
V2^, the part of the integrand that depends on s is just {zi/v2)'^. Therefore, 
deforming contours so that \v2\ > \zi\ gives ([5| and □ 



3.3 Analysis steps in proof of theorem 3.1 



For this section, we need a slightly different expression for the kernel. By (40)- 
(42) of [3j, the kernel equals 

K{ni,ai,si;n2,a2,S2,t) 



Je-'O J\z\ = l e*(^^T ' 

X ; q T^TT-a, (dz)dv 

(^+^' - 1)»2 z + z-1 - - -u-i ^ 



(31) 



+ I(ni,ai)t>(n2,a2) 



'y^(ai,-l/2)(g^) 



1^1 = 1 



-1/2) 



1(12,-1/2) 



-1 ma,{dz)\ (32) 



-1 \ ni-n2 
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'E^!iL!^!n£i)£"' j(..-v2) 



2 

ni — 712 



(dz) , (33) 



where 9 is any real number, and the arc from e to e*^ is outside the unit 
circle and does not cross (— oo,0]. 
Set 

z + z^^ f z + z~^ \ 
G{v,r],T,z) = T h?/log ( 1 j -vlogz 

By Proposition 5.1.1 of P], we can take V to be 

= {{'^,V,t) ■■■n,T> 0,(7i(77,t) <v < q2{ri,T)}, 

where 



qiir],T) = <^ 'V V V 2772 V T y ' - V "2 

1 T 




1-7:^ 1 + — , o<-< 



2-77' 

Lemma 3.3. Let fl± denote fl{±i/,ri,T). Then Ct^^l_ = 1. 
Proof. In general, 

r{z) 

where p and r are 

p(z) = r + (2?7 + 2i^ - t)z + (277 - 2i/ - r)z^ + TZ^, 
r{z) = 2z^{z-l). 

Let p±{z) denote the polynomial p(z) corresponding to (±7/, ?7,t). Note that 
z^p+(z~^) — p_(z). By definition, n± is the zero of p± that is in the upper 
half-plane. Therefore, flZ^ — □ 

Now let us return to the proof of the fourth condition in Definition |2.1[ 



Start by examining (31). Expanding the parantheses, we obtain four terms 
corresponding to z'^^v'^'^, z^^w"'*^, z"*!?;"'*^, and z~'^'^v''^. For the two terms with 
z'^'^ , make the substitution z ^ z~^ . What remains are two terms, corresponding 
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to z ^^v'''^ and z "^v Therefore, (31) equals 



l)"i 



1-v- 



dzdv 



l)n2 Z + Z' 



— V — V 1 2iz ' 



(34) 



We now need to deform the contours in (34) to steepest descent paths. In 
other words, we need 



(35) 

(36) 
(37) 

for all V on the w-contour. By Lemma [3. 3| and the definition of G, we see that 
?li{G{u2,m,r,fl{v2,rj2,T))) = 5R(G(~z.2, ?72, r, r!(-i/2, 7?2, r))). If \v\ > 1, then 
3?(G'(— 1^2, ?72, T, w)) > di{G{i'2,V2,T,v)). Since the steepest descent paths can 
go completely outside the unit circle (see Proposition 5.1.2 of [3]), (37 1 follows 
from (|36l). 



3?(G(i/i, 771, r,z)) < ^{G{iyi,r]i,T,Q{i^i,Tji,T))) 

for all z on the z-contour and 

3?(G'(z/2,'72,T,w)) > ^{G{iy2,V2,T,n{i^2,m,T))), 
n{G{-iy2,m,T,v)) > di{G{-i^2,V2,r,n{-i^2,V2,T))) 



Figure 5: On the left is 3? (G(i^i, 771, r, z) — G(i^i, r/i, r, ri(i/i, 771, r))), and on the 
right is 3?(G(i^2 7 v)—G{v2,il2,T, ^{'^2, V2,t))). White regions indicate 3? < 
and shaded regions indicate 3? > 0. The double zero occurs at Q{vj,rij,T). The 
arc V goes from e~*^ to e*^. The unit circle has been drawn on the right. 






n 






. 1 


{ 1 V 








\ 







If we deform the contours to the steepest descent paths Fi and F2 in Figure 
[5] we get that (31 1 asymptotically becomes 
2 



27ri 



cxp(jVG(t/i,7?i,T, z)) 1- v 
ri Jr2 CTq)iNG{v2,V2;T,v)) z + z-^ - v - v 



^ dvdz 

-1 -y 



22 



2ni 



exp(7VG(i/i,?7i,r, z)) 



1 - V- 



dvdz 



^ Jy'^ exp{NG{~V2,V2,T,v)) z + z ^-v-v ^ z 



plus possibly the residues at z = v. Since r2 goes outside the unit circle and 
the critical point of G(— 1^2, f?2, v) lies inside the unit circle, the second double 
integral is negligible. 

Now we need to compute the possible residues at z — v. If the contours pass 
through each other, then the residues at z = v equal 



^ , Z + Z ^ \ "1 »l2 ^ 



_1 X ni-n2 7 

z + z \ az 



(38) 



iy(-V2,-l/2)(s^) 



47ri 



dz 
z 



W^(-iA-i/2)(.,^) 



47ri 



-1 \ iii-ra2 , 

Z + Z \ dz 



, (39) 



where ^ is any complex number satisfying ( 35 1 and ( 36 1 . See Figure [6] If the 
contours do not pass through each other, then there is no contribution from the 
residues. For notational convenience, set 



C, if C exists, 
e'^, otherwise. 



Figure 6: The z and v contours from Figure [Sj They intersect at C,. 




It is important to note that ^ is arbitrarily selected. The only requirement 
on C, is that it satisfies the inequalities (35) and (36), and the only requirement 
on e*^ is that 3?(G2(e*^)) > 3?(G2(f^2))- So there exists e > such that if 
1^1 — ^1 < e, then ^1 also satisfies those inequalities. 

Now we need to compute (32) and (33). Expanding the parantheses, we 
get four terms corresponding to 2^1+*^ , z^^"^^ , z''^"''! , z"*"^"*^ . For the terms 
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corresponding z 



and 



*2 make the substitution 



Therefore, 



the sum of (l32|),(|33|),(138|),(139|) equals 



1 

47ri 



- 1 



1 

47ri 



dz 
z 

z" 



- 1 



dz 

z 



(40) 



where the contour crosses (0, oo) if n\ > n2, and it crosses (— oo,0) if rii < n2- 
For each integral, deform the contour to a circular arc of constant radius. Then 
the absolute value of the integrand is maximized at the endpoints. To sec this, 
note that \z\ is constant on the arc, so it suffices to analyze 



Rix,y) := 



{x + iy) + {x + iyY 



- 1 



Now note that 



Y7pr w A y{{x-l)^ + y^) 

Vi?(x, y) ■ (y, -x) = ^ . 

x^ + 



Therefore, R increases as we move counterclockwise along the arc in the upper 
half-plane, and R decreases as we move counterclockwise along the arc in the 
lower half-plane. If rii > ri2, then /j("i-"2)/2 maximized when R is maxi- 
mized, which occurs at the endpoints since the contour crosses (0,cxi). Like- 
wise, if rii < "-2: then R is minimized at the endpoints. Thus, in both cases, 
/j(ni-n2)/2 jg maximized at the endpoints. 

Using a standard asymptotic analysis (see e.g. chapter 3 of [ID]), we get that 
the asymptotic expansion of ( 40 1 is 



N 



- 1 



- 1 



- 1 



N 



N{r]i-r]2) 



N{rii-V2) 



- 1 



for some constants ci , C2 . To complete the proof, notice that if 



±1/2— 



- 1 



> e 



K(Gi(Oi)-G2(02)) 



for some selection of ±, then the asymptotic expansion of the kernel would 
depend on ^. But ^ was arbitrarily selected, so this is impossible. 

Now that the fourth condition has been proved, it remains to show that the 
second condition in Definition 2.1 holds. Recall that il{iy,ri,T) is the root of 
p(v, 77, r, z) that lies in the upper half-plane, where p is the polynomial from 
Lemma [3.31 We thus need to solve 



P('72(?7, t) - ei, 77, T, n{q2{'n, r), 77, t) + 63) = 0. 
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Since fl{q2{r], T),r], r) is a double zero of p{q2{'r], T),r], r, z), we thus have to solve 
^elp"{n) - 2ei(0 + 62 -el- 2e2fl - O^) + 0{el) = 0, 

1/2 

which implies that 62 = 0(ei )• In other words, as 1^ approahces q2{ri,T), 
Ct{u,ri,T) - Cl{q2{'n,T),r],T) = 0{{q2{r],T) - i^)^/^). Plugging this into the ex- 
pression for G" gives the result. 

4 Asymptotic Lemmas 

4.1 Riemannian Approximations 

Lemma 4.1. Suppose that g E [a, b] and I € [a, h] . Suppose thai as (5—^0, 
the Lesbesgue measure of the set {x € [a, b] : I'{x) € 27rZ + [—5, 5\} is 0{S°') for 
some positive a. Let e^ € [—1, 1] depend on N. Then 



[{b-a)N\ 

hm y e^^^(«+('=+^-)/%(a+^^) ^=0(1). 

fc=l ^ ^ 



Proof. Let <fe denote a + (/c + eM)/N. Note that — t^j = \k — s\lN. Fix some 
1 + A/'Vs < s < [(6 - a)iVJ - N'^/^ and consider 



fe = s-JVl/3 

We bound this sum in two cases. 

Case 1. Assume I'{ts) ^ 27rZ + [-5,6]. For s - N^/^ < k < s + N^/^, 
Taylor's theorem says that 

I{tk) = [I{ts) + I'{ts){tk - ts)] + [ll"{ck){tk - tsf] =: [Ii{tk)] + Mtk)] 
for some Cfc between tg and tk. We will prove that 

2^ e g [tk) ^ S ^ + ^ + ^4/3 

k=s-m/3 



Using the inequality 

\9{tk) - g{ts)\ < \\g'\\oo-\tk-ts\ = y 

we have that 

< 

N 



I \k-s\ 

1 00 



N 



/s=S-JVl/3 



< E ll^'lloo^ 

fe=S-JVl/3 
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Furthermore, for \k — s\ < 7V^/^, 

|l_e«JV/.(t.)| ^ |i_g./"(c.)(fe-^)V(2JV)| < <9||/"||^Ar-i/3 
Also, 



(42) 



A;=s-Ari/3 



„j/'(t,)(/c-s) 



< 



< 995-^ (43) 



Using (|41|, the definitfon of Ii and /2, (|42|) and (|43p respectively, 
1 



E 



TV 



< 



s+N 



1/3 



'(^«) E 



9[ts) ? e 



,1 

N 



+ 3||5'||ooiV-'/' 



fc=s-Arl/3 



E 



fc=s-Arl/3 



N 

18||.g||oo||/"||oo , 3||g' 



3||ff'||oc 



N 



< 



99S-^g\\^ , 18||5lloo||/"||cc. , 3||g'||, 



N 



N 



Case 2. Assume that I'{ts) E 2ttZ + {—S,S). In this case, only a simple 
estimate is needed: 



fc=s-Ari/3 



< 



Ar2/3 



Since the estimate in case 1 is much better than the estimate in case 2, we 
need an upper bound on how frequently case 2 can occur. In other words, we 
need an upper bound on the measure of the set {a; € [a, b] : I'{x) G 27rZ + (5, d)}. 
We assumed that 

\{x G [a, b] : I'{x) G 27rZ + [,5, = OiS"-). 

Now we need to sum over aU s in the set {N^/^ + l,3iV^/^ + 1,57V^/^ + 
1 . . . , (6 — a)N — iV^/'^}. There are 0{6''N'^/^) terms for which case 2 applies. 
Therefore, 



L(b-a)ArJ 

E 

k=l 



JNI(a+{k+tM)/N), 



k + 1 



TV N 



and setting S = N yields the result. 



□ 
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Proposition 4.2. Suppose f : Z>o — > C is a function such that for each t > 0, 
f{ltN\ ) = e''^''^*^g{t)N'^ + o{N'^) as TV ^ oo, 



where g and I satisfy the same assumptions as in Lemma \4-l\ Further suppose 
that the error term o{N'^) is uniform, i.e. 



f(ltN\)-e''^'^'^g{t)N'^ 



— ^ uniformly on [a,b]. 



Then as N ^ oo, 

mi 

x=[aN]+l 

Proof. Let e > 0. By assumption, there exists Nq such that for N > Nq and 
xe [laN\ +l,lbN\] 

m - e»^^(-/^)g (^) N^^ 

By summing over x, we see that 

ybN\ 

x=[_aN\+l 

The Proposition now follows from Lemma |4.1| 



d+l 



< {b~a)eN 



□ 

Lemma 4.3. Suppose g is a function on [a, b] with bounded variation. Then as 
N ^ oo, 



E 3{^)=N g{t)dt + o{N). 



x=[aN]+l 

Proof. Let e > 0. By the definition of a Riemann sum, there exists M such that 

L(6-a)7VJ 



/ 9{t)dt- E 
k=i 



g[a + 



k\ 1 



N N 



< e/3 for aU N > M. 



Furthermore, for N > 3BV{g)/e 

[ib-a)Ni 



- E 



5 a + 



k {aN} 



Finally, for N > 3||5||oo/e 

[{b~a)Ni 

N 



- y 



k=l 



k \aN} 

g I a H -'^ - 



[bN] 



- E 



< 



□ 
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Proposition 4.4. Suppose f : Z>o — > M is a function such that for each t > 0, 

filtN\ )=g{t)N'^ + oiN"^) as N ^oo, 

where g is a function on [a, b] of bounded variation. Further suppose that the 
error term o{N'^) is uniform, i.e. 

f(\tN\) - g(t)N'^ 

> uniformly on [a,b\. 

Then 

J2 fix) = N^^^ I g{t)dt + o{N'^+^). 

x=[_aN\+l "^^ 

Proof. Let e > 0. By assumption, there exists A^o such that for N > Nq and 
xe [laN\ +l,lbN\] 

g (I) iV^ - eN'' < fix) < g (|) iV^ + eN''. 

We can see this by noting that when t = {[aN\ + 1)/N, [tN\ = [aN\ + 1. 
Furthermore, increasing x by 1 causes t to increase by which causes ltN\ 
to increase by 1. By summing over x, we see that 

lbN\ lbN\ lbN\ 

2;=LaA'J+l a;=LaArJ+l x=LaJVJ+l 

The Proposition now follows from Lemma |4.3| □ 



4.2 Asymptotics 

Proposition 4.5. For j = 1,2, let {vj^rjj^r) E T), Vlj denote Vl{vj,rij,T), 
Gj{z) denote G{vj^rjj,T,z), and 0j denote 6{i>j,rij,T). With the assumptions 
from section\2.1\ 



l \ f f exp{NG{rii,vi,T,u)) 

' ' ' f(u,w)dwau 



, 3 -L n"fa.\-3 

O 



G'd^i) + 02(^2) j^-2\ ,n(r"(ri \-T/'^r"(rt ^-V2A^-3^ 

„ArsR((Gi(Oi)-G2(02))) r iAfa(Gi(Oi))-i9i _ iAra(Gi(ni))-ifl 



2nN^\G'i{n,)W\G'i{n2)\ 



iiV3(G2(02))+ie2 ■' ^ g-iA'S(G2(02))-ie2 



-iArs(Gi(ai))+iei _ _ p-jArs(Gi(ni))+iei 



giArS(G2(02))+ie2 •'^ ^^g-iAr9'(G2(n2)-i6l2) 

Proof. First, we show that the main term is correct. 

By assumption, we can deform Fi and F2 so that Fj passes through flj, Qj for 
j = 1,2. The contributions to the integral away from flj^rij are exponentially 
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small, so we can replace Tj with 7^ U jj, where jj and are steepest descent 
paths near fi^ and Clj, respectively. The integration over u G 71U71, u; G 72U72 
expands into four integrations corresponding to {u,w) £ 71 x 72,71 x 72,71 x 
72, 7i X 72- We explicitly do the calculation for 71 x 72. The other three 
calculations are essentially identical. 

Make the substitutions s = Gi{fli) — Gi{u) and t — 02(^2) — G2{w). In 
the neighborhood of m = 57i and w — Vl2, we have 

/(u,«;)«/(Oi,172), s = -^^^^^G'/(l^i), i = -(^^^^^G'2'(l)2), 
which imply 

G'M = - J = - ^i)G'i{^i) = ^-2sG'i{n,), 
dt 
dw 

Then we get 



G'2{w) = -^ = {w- n2)G'^{n2) - ^-2tG'i{Sl2). 



Af(Gi(ni)-G2(02)) ^00 .00 

4- : ^ f{n,,n2)[j s-'/^e-^^ds)[j r^'^-'^'dt) 



87r^y^GU^)y^GWh) 



gAf(Gi(ni)-G2(n2)) 



giV3?((Gi(Oi)-G2(02))) 



27riVV|G'/(f7i)|V|G^'(^^2)| 



/(^^l,f^2 



giAfa(Gi(Oi))-iei 



giAfS(G2(02))+ie2 



where the last equality follows from G{z) ~ G{z). The 4 appears because the 
maps M I— ^ s and w ^ t arc both two-to-one. 

It still remains to show that the error term is correct. The remainder of this 
section is devoted to proving this. The idea is to reduce the double integral 
to progressively simpler forms. First, by a reparametrization, the integral over 
two arcs in C can be written as a integral in M^. Second, by using a Taylor 
approximation, the integral in can be written as a product of two integrals 
in M, each of which is of the form / e~^^'-*-'0(i)di, where R{t) has a maximum 
imax in the interval of integration. Third, by using the implicit function theorem, 
this integral reduces to the form j e^^* where the interval of integration is 
a small neighbourhood imax- Fourth, this last integral is a slight generalization 
of / e~^^ g{t)dt, which is dealt with by the well-known Watson's lemma (Lemma 
|4.6| below). Since the first two steps have been done before (see Chapters 3 and 
4 of [TU]), we will focus mostly on the third and fourth steps. 

Lemma 4.6. Suppose T > and 4>{t) is a complex valued, absolutely integrahle 
Junction on [0, T] : 

rT 

\<l){t)\dt < 00. 
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Also suppose that (j){t) = g{t), where cr > — 1 and g is continuously dijferen- 
tiable in some neighbourhood oft^O. Then for any N > 1 and s S [0,r], 



r(<7 + 2) sup \g'{T)\ 



m)\dt 



+ 5(0) / e-^H^dt + e~^'/^r{3 + 2a) sup |.g'(r) 

0<T<S 



Proof. See Proposition 2.1 of [lOj . 



□ 



Lemma 4.7. Suppose T > and 4i{t) is a complex valued, absolutely integrahle 
function on [0, T] : 

\(l){t)\dt < oo. 

Also suppose that (f> is twice continuously differentiable in some neighbourhood 
oft — 0. Then for any N > 1 and s e [0, rri^], 



-Nt' 



cj){t)dt 



iVl/2 



< 



sup \g'{T)\ 

'TT 0<r<s 
2 WJ^ 



m)\dt 



\(f>{t)\dt + (l,{Q) I e-^H-^'^dt + e'^''^ sup |.g'(T)|, 

0<r<s 



where m = inin(a,/3) and 
Proof See pages 58-60 of HUj. 



□ 



Lemma 4.8. Suppose that R and (p are infinitely continuously differentiable in 
some neighbourhood of tmax- Also suppose that tmax *s a local maximum of R 
and R"{t,nax) < 0. Then for any N > 1 and s G [0,m^], 



+ e 



where 



tmoi-<5l 

Ns 



^ max(Q,/3) 



-2n 



NR"{traax) 



< 



sup |g'(r) 
V^r o<T<s 

~2 7V372 



\h{t)\dt + e 



\h{t)\dt + h{id) / e-^H-^'^dt 

J S 

sup |5'(r)|, 



0<r<s 



= \/-R{tmax - h), P = \/~R{t,nax + h) , m = mm{a, /3), 
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his) = + sv{s)){sv'{s) + v{s)), g{s) = ^{h{s^/^) + h{-s^/^)), 

where v{s) is an infinitely differentiahle function solving 

~R{tmax) + R{tmax + Sv{s)) = - S^ . (44) 

Before continuing, a few estimates on v{s) are needed. 



Lemma 4.9. Let v{s) be as in p4[ ). 

(a) 

R"{t,„ax) = -2vi0)-^. (45) 

v'iO) = ^"'^^^--\ (o)^ (46) 
(b) Set B = sup |i?W/24|. Then 



_ ^,(0) - v'{Q)s\ < |^-^i?"'(t„,,)2„(0)7 + BviOf ] s\ 



for 

s < mm{\53R"'{t„,ax)\/{625Bv{Q)), {\R"' {t^a.)HOf)-\ {50VBv{0f)-^). 

In particular, \v{s) - v{Q)\ < \R"' {tmax)\v{Oy\s\/'l and \v{s) - v{0)\ < u(0)/4. 
(c) Let 



Then 



471 

\v{s) + sv\s) - v{0) ~ 2v\0)s\ < { gi?"'(t,„,j2„(0)7 + -^BviOY 



s 



for \s\ < min 



|53i?"'(t. 



'2w(0)-i 



3a3 



625Bt;(0) ' 50/Bt.(0)2 ' 6i?"'(i™,,)«(0)= 

' (47) 

(d) With the same bounds on \s\, 

\2v'{s) + sv"{s) ~ 2w'(0)| < 450000(i?"'(t™a:,)2w(0)^ + Bv{Gf )s 
Proof, (a) The proof comes from page 69 of [10 . It follows immediately from 



using implicit differentiation of (44) and setting s — 0. 

(b) First notice that if i?_(t) < R{t) < R+{t) with i?_(imax) = ^(^max) = 
i?+(tniax) and v± are the solutions to — i?(<,nax) +^(imax ± si'±(s)) — — s^, then 
V- < V < v+. We will use 

R±{t) = i?(imax) + 2^"(^max)(t — imax)^ + g (^max) — ^max)^ ± B {t — t„ji^y^)'^ 
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Therefore, wc obtain bounds on v{s) by solving — ii(fmax) + -R(tmax±si'±(s)) = 
— s^, which is equivalent to solving 

<3±,«(2/):=l-2/oV + ^S2/3±BsV = 0, A = R"'{t^^)/6, yo=viO). 

In other words Q±,s{v{s)) = 0. We will use the intermediate value theorem to 
estimate roots of Q±^s- 

For e = Ay^s/2 + {bA^yl/A + ByDs^ and |s| = {HAyl)-'^ where H is any 
real number, we have 

Q+AV^ + e) < 256AiogiOyio (256g^^ + 2bQA^ B^p^{H)yl + 32^1453^4 (h)^^ 

where the Pi are polynomials which satisfy the following inequalities when | > 
6 

P2{H) = {5 + 2H + 4H^) < IIH^ 

Pi{H) = (75 + mn + 132i/2 + 56i7'^ + A^H'^) < 371H^ 

Pq{H) = (125 + 150ff + 360i?2 + 308i/^ + 312i?* + lUH^ + 48H^) < 1U7H^ 
Ps{H) = (625 + lOOOiJ + 2600i/2 + STGOH^ + 4336i7* + AIQOH'' + 1792il6 

+ 1024il^ - 256i?^) < -lOH^ 
H^P5{H) = H^{125 + 150H + 360H^ + 148H^ + 208H^ - 80H^) < 

Now setting H := hyQ^A^^^B where h > 50, 

^ , , 128 + 1408/i2 + 5936/1-* + 11576/i6-5/i« „ 
Q+AVO + e)< < 0. 

Since 

Q+Ayo + Ayts/2) > ^s^y^{B{2 + Asy^o)^ + 2AYo{lO + 6Asyl+A^s^y^o)) > 0, 

this implies that v+{s) < ya + e = v{0) + v'{0)s + {5A'^v{0y/4 + Bv{Of)s'^ for 
|s| < min((6mt;(0)-'')-i, (50/Bt;(0)2)-i). 

By an identical argument, for e = AyQs/2 — (5A^j/q/4 + i3j/o)s^ and \s\ = 
{HAy^y^ where \H\ > 6, 

Q-Avo+e) > 256^iOffio^io (-256i?^^256A^i3^(llg^)yg-32^^i?^(347g^)^/o" 

- 16A^B2(1167iJ^)i/^ + A^B{86H^)y^), 

and setting H := hyQ^A~^y/B where h > 15, 

^ , , -128 - 1408/i2- 5552/i4-9336/i'5 +43/1^ „ 
Q-Avo + e)> > 0. 
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For As>0 and \s\ < niin(|318A|/(625Syo), {6\A\y^)-^), 

g_,,(2/o + Ay^s) = si-AyU2 + Asy^) + Ay^^i^ + Asy^f - Bsy^{l + Asy^)^) 

s{318A + 625Byos)y^ 
1296 

Therefore yo + Ay^s/2 - {^A^yl/'i + ByDs^ < v^{s) < yo + Ay^s for As > 0. 
For As <0 and \s\ < \A\/{Byo), 

Q-Avo) = syliA- Bsyo) < 0, 

so 2/0 + ^2/0^/2- {5A'^yl/4: + By^)s^ < v^{s) < yo. Thus the lower bound holds 
in both cases. 

The last statement follows because 

Hs) - vm < \v'm ■ \s\ + (^A-R"'(t^^,yv(oy + BviOf'^ s' 

< 12 l^"+144|i?'"(W)K0)3"^"+^"(°) 625M0) 

^ |i^"'(W)b(0)^ , , ,^(0) 
4 4 ' 



(c) Differentiating (44) yields 

-2s 



or equivalently 



Hs) + sv'{s)) = - 



N -2s - v{s)R'{tm^^ + svjs)) 



To estimate v'{s), let us first estimate R' . 
By a Taylor expansion, 

|i?'(imax + SV{S)) - " ^ i?'" (tmax)s'«(s)' | < 4Bs3„(s)3. (49) 

By the triangle inequality and part (b), 

i?'(imax + SV{S)) - i?"(t„,ax)s(«(0) + Sv' {0)) - ^R'" it,,.^)s\{0)^ 

< ABvisfs"" - - viO) - v'{0)s)s + - v{0)){vis) + v{0))s^ 

< ^BviOfs^ - R"{t^,^) (^±^R"'{t^,^)^v{Oy + BviO)''^ s' 

+ Ji?"'(W)W•^«(0)^^ 
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which, by (45) and (46), implies 



|i?'(imax + sv{s)) + 2vi0)-\s - 4t; (0) "^Z (0) s 



< 



16 ' ' 288 
To estimate the inverse of R' , use 



--■.a3s\ (50) 



1 



1 



a2 



{al - aia3)s^ + 02033^ 



al{ais + 025^ + 035^) 
for \s\ < min( 



< 



6(a| + laiasl) 



ais + a2S^ + a^s'^ ais 



which, by setting ai = 2v{0) ^ and 02 = 4i;(0) ^w'(O), imphes that 
1 i;(0) 







3a2 





R'{t 

max 

+ sv{s)) 2s 



v'{0) 
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Multiplying by 2|s| finishes the proof of (c). 
(d) Differentiating ( 44 ) twice yields 



2v' (s) + sv" {s) 



-2 - R"{t^^^ + sv{s)){v{s) + sv'{s)f 



(52) 



-R'(imax + SV{S)) 

From part (c) and a Taylor approximation for i?", 

I -2 - i?"(W + sv{s)){v{s) + sv'{s)f\ < 999((i?"'(W)'i^(0)^ + BviO)')s' 
+ {R"'{U..fv{Of + R"\t^,^)Bvi0Y)s' 
+ (i?"'(W)"«(0)i2 + i?"'(Wx)'i?«(0)i° + mB'v{Of)s^ 
+ (i?"'(W)'f (0)^^ + i?"'(Wx)'St'(0)^^ + 10i?"'(W)B'«(0)")s5 
+ (i?'''(W)^(0)i«+i?'''(W)'i?«(0)i^+10i?'''(W)'i3'«(0)i^+10i?3^;(0)i2)s6). 

Since |s| < {R"'{UMOfr\R"'{tm..){Bv{0)r\ 

\-2 - i?"(Wax + sv{s)){v{s) + sv'{s)f\ < 99999(i?"' ( W)'«(0)' + BviOf)s^. 

(53) 



By (51) and the estimates on |s|, 

1 v{0) 



max 

+ sw(s)) 2s 



«'(0) 



< 30i?"'(tn,ax)«(0)^ 



(54) 



Combining (l52j),(l53|) and (154|, 

\2v'{s) + sv"{s) - 2u'(0)| < 50000(i?"'(tmax)^w(0)^ + Bv{Qf)s 

+ 40oooo(i?"'(w)^«(o)i° + i?"'(^,„ax)B«(o)«)s^ 

and using \s\ < (i?"'(tmax)w(0)'^)^^ on the second term gives the result. 



□ 
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Corollary 4.10. Suppose that R and (j) are infinitely continuously differentiable 
in some neighbourhood oftmax- Also suppose that t^ax is a local maximum of 
R and R"{tmax) < 0. Let 6i and 62 be positive numbers such that 

-RiUnax - Si) = -R{tmax + 52): 

and assume equals the right-hand side of (|47|) . Let 



L^ ij^max) 1 
50Bw(0) ' bm"'{tmax)v{QY ) ' 

A := 5OOi?'''(i„,J||0'|Ut'(O)'V45OOOO||(/.||oo(i?'"(tmaJ'i'(O)' + Sz;(0)5). 
Then for any N > 1, 



tmoi-<5l 



/ -27r 

NR"{t^ax) 



< 



Proof. Use Lemma 4.8 By part (a) of Lemma 4.9 



By parts (c) and (d) of Lemma 4.9 



sup |ff'(T)|<A. 

0<T<m2 

When t;(0) > 

53 i?"'(i„,ax) ^ 53 



> 



625 Bv{0) 625i?"'(tniax)w(0) 

1 1 

> 



3 



50/B«(0)2 50i?"'(W)w(0)3 
as < l&R"'{t^^^fv{Qf 



/2w(0)-i ^ 1 



y 3a3 5i?"'(W)i'(0)3- 
When v{Q) < 

1 ^ -R"'(imax) 



i?"'(imax)^'(0)3 S«(0) 
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1 ^ -R"'(imax) 



as < l6Bv{0f 



Thus > s, and 



3a3 5/Bi;(0)2 5Bv{0) 



mN y/SN' 

□ 



We can finally wrap up the proof of Proposition [15] Since s is not too small 
(polynomial in v{0)~^), the exponential terms are small enough to be ignored. 
Therefore the error term is Ai = O{v{0y N"^/^) = 0{R'{{tn,i,,,y^^^N'^^^). 
The main term is of order N~^/^R'({tma.yi)~^^'^ ■ Thus, when multiplying two 
integrals of the form in Corollary |4.10[ we get 

as needed. □ 



In Proposition |4.5[ the error term blows up at the edge. Therefore a better 
bound is needed. To get this bound, we simply use the first term in Watson's 
lemma, as opposed to using two terms. Since the method of the proof is identical 
as before and the details are simpler, the proof will be omitted. The exact 
statement is the following. 

Proposition 4.11. For j = 1,2, let (i/j,rij,T) E T), Vlj denote Vl{vj,rij,T), 
Gj{z) denote G{i>j,rij,T, z), andOj denote 9 (vjyrjjjT) . With the assumptions in 
section \2.1\ 



l \ f f exp{NG{r]i,iyi,T,u)) 

' ' ' j(u,w)dwdu 



2mJ Jn Jra exp(iVG'(7y2, 1^2,1", 
1000 



|/(i]i,r!2)|+|/(r!i,n2)|+|/(JTi,r!2)|+|/(JTi, 172)1 



A Gaussian Free Field 

This section is adapted from jllj . 

Given a domain D C M'*, let Hs{D) be the space of smooth, real- valued 
functions on M'' that are supported on a compact subset of D. Let (f,g)v = 
/^(V/ • Vg)dx be the Dirichlet inner product on Hs{D). Let H{D) be the 
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completion of Hs(D) under the Dirichlet inner product. Note that H{D) is a 
Hilbert space. 

A Gaussian free field is any Gaussian Hilbert space Q{D) of random 
variables denoted by {h, /)v - one random variable for each / e H{D) - such 
that 

E[ih,a)v{h,b)v] = {a,b)v- 
In other words, the map from H{D) to Q{D) defined by sending / to the random 
variable (/i, /)v preserves the inner product. 

By the identity (a, 6) = ^[{a + b,a + b) — {a,a) — {b,b)], this map preserves the 
inner product if it preserves the norm. In other words, the variance of {h, /)v 
is (/,/)v for each f€H{D). 

At the most fundamental level, a Gaussian free field is a Gaussian Hilbert 
space. It turns out to be equivalent to consider other Hilbert spaces, which we 
will consider below. 

Let {cj} be the eigenfunctions of the Dirichlet Laplacian on D which form 
an orthonormal basis for L'^{D), and let {Xj} denote the corresponding eigen- 
values. By the spectral theorem for compact self-adjoint operators, all the Xj 

are negative. Suppose the ordering satisfies Ai > A2 > Define the operator 

on L'^{D): 

(-Ar5^/3,.e,.=^(-A,.r/3,.e,.. 

Define Ca{D) {~A)''L'^{D) to be the set of all ft e^ such that E(--^j)""/3iej € 
L'^{D). We may view Ca{D) has a Hilbert space with inner product {f,g)a = 
((—A)""/, (— A)""^). Note that when a < 0, the sum ^ pjCj may not converge 
in L'^(D), but it converges in the space of distributions. 

Note that by integration by parts, H{D) = £_i/2{D). Therefore, instead of 
considering H{D) as the Hilbert space, we can also consider £1/2 = (~^)^(^)- 
Then 

E[(/j,/)i/2(/i, 5)1/2] = (/, 5)1/2 = {{-A)-^/'f,{-A)-^/'g) = ((-A)-V,5). 
Since 

[(-A)-V](x)= / G{x,y)f{y)dy, 



D 



this implies 

Cov[(/i,/)i/2,(/i,5)i/2] = / f{y)g{x)G{x,y)dxdy. 

JdxD 

The higher moments can be computed from the following theorem. 

Theorem A.l. IfXi, . . . .X^ are mean zero random variables such that aiXi + 
. . . + akXk is Gaussian for all ai, . . . ,ak &M., then 



k/2 

^ J|e[X<^(j)X^(j-+i)], k even 
0, k odd, 



nxi...Xk] 

where the sum is over all involutions on {1, . . . , fc} with no fixed point. 
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